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Abstract 

The strength function of a single state \d) is studied using the deformed 
Gaussian orthogonal ensemble. In particular we study the dependence of the 
spreading width of | d) on the degree of mixing. 



The mixing of a single state with a background of complicated states is important in the 
description of a variety of phenomena in nuclear physics such as isobaric analogue resonances, 
giant dipole resonances and the decay out superdeformed rotational bands. Such can mixing 
be conveniently described by the strength function [0] and it is interesting to study the 
generic features of this object. Previous studies in this vein have investigated a single state 
coupled to a background generated by a two-dimensional anharmonic oscillator p| and the 
spreading of a shell model basis state over the shell model eigenstates due to the the residual 
interaction 0. Here we use random matrix theory [Q. 

We write the Hamiltonian, H, as the sum of two terms, 

H = Ho + V. (1) 
The eigenstates, |n), and eigenvalues. En, of H satisfy 

H\n) = En\n), n = l,...N + l, (2) 



whilst for Hq we have 



Ho\k) = Ek\k), k = l,...N, 

Hold) = Ed\d), d= N+1. (3) 



The strength function is then defined as 



Af+l 

F,{E)=Y: \{d\n)\'5{E-En), (4) 

n=l 



* Supported in part by FAPESP. 
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and describes how the state \d) is distributed over the + 1 eigenstates of H. We shall 
choose the energy of \d), E^, such that it lies in the middle of the the spectrum of H. 

Some insight into the behaviour of the strength function can be obtained by performing 
a two step diagonalisation of i^f Let us represent the Hamiltonian in the basis \d)}. 
Diagonalising H in the A^-dimensional subspace defined by excluding \d) we obtain the set 
of eigenvectors \q) = J2k{k\(l)\^) with eigenvalues Eg, q = 1, ...,N. The Hamiltonian in the 
basis {l^z), has diagonal matrix elements Eg, q = 1, N and Ed, d = N + 1 and non- 
zero off-diagonal elements Vdg = Vgd = J2k^dk{k\Q) (the Hamiltonian is assumed to be real 
symmetric). The diagonalisation of the intermediate matrix may be carried out analytically 
|Jl|] so that using a Lorentzian of width / to represent the ^-function in Eq. (^), one obtains 

i 



where 



and 



(E - E,f + 



2' 



By making the further assumptions that the eigenvalues Eg are equi-distant with mean 
spacing D, that the squared matrix elements iVg^p have approximately the same order of 



magnitude, (V^), for all q and that the magnitude of y (V^) is smaller than the energy range 
in which it may be considered constant (whilst being larger than D in order for the strength 
function to have meaning) the strength function may be approximated by a Lorentzian 

^'^(^) ^ ^ p. . (8) 

2^ {E - Edf + i^-if 



)5 



where the spreading width is given by the "golden rule 

ri-2vri^. (9) 

We wish to study how the distribution, Fd{E), depends on the degree of mixing by 
which we mean the strength of V. To that end we employ deformed Gaussian orthogonal 
ensemble (DGOE) In this model H is real symmetric and it's matrix elements are taken 
to be independent Gaussian distributed random numbers with zero mean and variances 
((//o)fcfc) = ^ for the diagonal matrix elements and ((V)^, = for the off-diagonal 
matrix elements. We take ((V)^^) = ((V)^^) = as well although choosing a different 
variance for this matrix element may be appropriate in some applications. 
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The parameter A may be varied between and 1 and determines the degree of mixing. 
The parameter a determines the energy interval over which the eigenvalues are distributed. 
In the limiting case A = the eigenvalues are Gaussian distributed; 



whilst in the limit A = 1 they are distributed according to the Wigner semicircular law; 

pi(^) = iV^Va2 - E-^. (11) 

We perform an unfolding of the spectrum of the spectrum of H defined by 

x„= [^"dEp{E), (12) 

where p{E) is the smoothly varying part of the level density. Thus the unfolded spectrum 
has mean level density equal to 1 and is dimensionless. The smooth variation the level 
density was obtained in practice by fitting the cumulative level density (staircase function) 

N+l 

x=Y.e{E- E„) (13) 

n=l 

(the letter 9 denotes the unit step function) to a polynomial using the method of linear least 
squares 0. For the case N = 50 (used in the calculations below) we found the function 
which gave the best visual fit was a polynomial of degree 5. 

In order to compare our calculation of the strength function on the unfolded energy scale 
with the Lorentzian approximation, Eqs. (^ and (^, we apply a simplified unfolding to 
these approximate formulae, defined by 

X = |. (14) 
Thus the unfolded version of the golden rule is 

D " (^^^ 

Persson and Aberg proposed the following formula for the mean (meaning averaged 
over the energy range where the eigenvalues are distributed) density of states of the DGOE 
which interpolates between A = and A = 1, valid when a = 2: 

PA = 1 . (16) 
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FIG. 1. Strength function, F^^E), for various values of the mixing parameter, A, calculated 
using DGOE (solid lines) and calculated using the Lorentzian approximation (dotted lines). See 
text for discussion. 



Figure |T] displays our calcuations of the strength function (solid lines) for various A using 
Eqs. (1^-0) where the eigenvalues (unfolded using the procedure defined by Eq. ([12|) ) and 
eigenvectors are generated by the DGOE. We performed our calculations using = 50 
whilst the width of the the Lorentzian used to represent the (5-function (Eq. (D) is / = 3. 
We set a = 2. 

The dotted lines in fig. |l| are the strength function calculated using the Lorentzian 
approximation, Eqs. (|^) and (^). The mean square value of the coupling is taken to be 
(V^) ~ The mean level spacing D is taken to be ^^^^^^^ = 2-^; the lower hmit for the 
level spacing of the DGOE spectra. 

Using the estimates of the previous paragraph, the condition > 1 (for the single 
state \d) to be significantly mixed with more than a single one of the \q)) implies A > 
Thus for = 50, although we can calculate the strength function for arbitrary small A it 
is only meaningful for A > 0.07. The value of / should be chosen so that it is a negligable 
fraction of the combined width F;!; -|- / whilst being greater than the level spacing (unity 
after unfolding). 

From fig. |l] we can see that for weak A the strength function has an approximately 
Lorentzian shape. As we decrease A below 0.07 the strength function approaches a Lorentzian 
whose width is dominated by I. For A > 0.07 the strength function starts to deviate 
significantly from the Lorentzian shape already at A = 0.1 broadening towards a semicircle 
before A = 1. An ensemble average was performed over 100 realizations for the cases 
A = 0.05, 0.08 and 0.1 whilst to obtain a relatively smooth strength function for the A = 0.2, 
0.5 and it was necessary average over 1000 realizations. 
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FIG. 2. Spreading width of \d) as a function of the mixing parameter A. See text for meaning 
of (iHiv). 

In fig. ^ we plot the following calculations for the spreading width of \d) as a function 
of A: 

(i) Eq. (H), at the peak value E = Ed where the eigenvalues (unfolded using the procedure 
defined by Eq. ([T^)) and eigenvectors are generated by the DGOE. An ensemble average 
was performed over 100 realisations for all A; otherwise the same parameters are used as 
were used in fig. |l]. 

(ii) The limiting value (A = 1) for the width calculated using the golden rule (Eq. |15|) and 
density Eq. (O) at i? = 0: 27r (^) (^)' = ^. 

(iii) The golden rule expression for the width using the Poisson density, Eq. (p!OD, at i^^ = 0: 

')nT(>^\ ( - ^ 

(iv) The golden rule with the density Eq. (|T^): 2ti p\- 

(v) The full width at half maximum FWHM of the strength function Ffi{E), Eqs. (§|-0), 
calculated using the DGOE. 

Three regions may be identified in the DGOE calculations (lines (i) and (v)) in fig. 
We see that T^diEd) for very weak coupling (A < 0.04) has a quadratic dependence which 
accurately follows the golden rule calculation (line (ii)). Between A = 0.05 and A = 0.15 the 
this dependence is linear becoming very weak as A approaches unity. The behaviour of the 
DGOE is simulated by Eq. (|16D of Persson and Aberg for the average level density (line (iv)). 

At A = 1 our DGOE calculation of Tli{Ed) is approximately a factor (^^^^^^^ = (^^^ = 1.6 
greater than the calculation which employs Eq. ([T6| ) . The FWHM (line (v)) also has a 
quadratic dependence on A for weak A. A linear dependence is maintained up to A = 0.3 
after which the A-dependence becomes very weak as A approaches unity. When A = the 
FWHM is just equal to the averaging interval /(= 3). In the opposite limit of A = 1 the 
FWHM is essentially equal to A^(= 50); ie. the state \d) is spread over all the eigenstates of 
H. Thus in the limit A = 1, the golden rule (line (v)) is a factor ^ smaller than the FWHM. 
Ref. 1^ investigated the spreading width of a basis state of a shell model Hamiltonian 
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as a function of the strength of the residual interaction (corresponding to Xa above). They 
identified regions where the spreading width has a quadratic dependence on the A (weak 
mixing) becoming hnear for larger A. They also found that the golden rule cannot be used 
to estimate the FWHM of the strength function for strong mixing. 

In conclusion we have studied how the shape and width of the the strength function of a 
single state depends on the degree of mixing, using random matrix theory. An application 
of these results to the decay out of a superdeformed rotational band is in progress. 

A.J.S. acknowledges helpful comments on the manuscript by Professor R.C. Johnson. 
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